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Abstract 

A variationally improved Sturmian approximation for solving time-independent Schrodinger equation 
is developed. This approximation is used to obtain the energy levels of a quartic anharmonic oscillator, a 
quartic potential, and a Gaussian potential. The results are compared with those of the perturbation theory, 
the WKB approximation, and the accurate numerical values. 
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' Since the early days of quantum mechanics, the main technical tools for solving the time-independent Schrodinger 
' equation have been the time-independent perturbation theory, the semiclassical or WKB approximation, and 
the variational method j^, ||. Starting form the late 1950's, physical chemists and nuclear physicists have ex- 



pH . plored the use of what is called the Sturmian basis functions in solving this equation for a variety of potentials 
arising in molecular and atomic physics ^. Recently, Antonsen j^] and Szmytkowski and Zywicka-Mozeiko 
[H have studied the harmonic oscillator Sturmian functions. The purpose of the present article is to outline 
a general variationally improved Sturmian approximation scheme that provides a nonperturbative method of 
, solving time-independent Schrodinger equation. 
^ , The organization of the article is as follows. In section II, we give the definition of the Sturmian basis vectors, 

derive their general properties, and discuss the conventional Sturmian approximation. In section III, we present 
an improved Sturmian approximation which makes use of the variational method. In section IV, we study the 
harmonic oscillator Sturmians and use them for the solution of time-independent Schrodinger equation in one 
dimension. In section V, we apply our general results to some concrete problems, and compare our results 
with those obtained using perturbation theory, the WKB approximation, and the highly accurate numerical 
investigations. In particular, we obtain the energy levels of a quartic anharmonic oscillator, a quartic potential, 
and a Gaussian potential. Finally, in section VI, we summarize our results and present our conclusions. 

II. Conventional Sturmian Approximation 

Consider the time-independent Schrodinger equation: 

H\E,a)^E\E,a) , (1) 
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where H is a self-adjoint Hamiltonian operator and a is a degeneracy label. 

The method of Sturmian approximation is based on an expansion of the eigenvectors \E,a) in terms of 
solutions \4>u,ct) of the equation 

>, (2) 

where Hq and Vq are self-adjoint operators, /3y and £ are real scalar parameters, and a is a degeneracy label. 
Note that in order to obtain a), one must fix £ and solve Eq. for \(j)^,a). Clearly, every solution a) 
would correspond to a choice for the value of the 'coupling constant' /3i,. 

Suppose that Hq = jfi /{2'm) is the Hamiltonian for a free particle moving in the configuration space H'', 
Vo = Vo(af) is a real interaction potential, and iJ is a standard Hamiltonian of the form 

i7 = |^ + y(f), (3) 

where p and x are momentum and position operators, respectively. If Vo(a^) tends to infinity as |a;| — > oo, all the 
eigenstates of Vq are bound states 0. In this case, only for a discrete set of positive values of can we find 
square-integrable solutions \4>v,oi) of Eq. (Q). In this case, the label v will take values in a discrete set which 
we shall choose to be {0, 1, 2, • • •}. 

The only difference between Eq. (|^) and the eigenvalue equation for the potential l3yV{){x) is that in the 
former f is a fixed parameter which can be arbitrarily chosen. Therefore, a solution of Eq. (H) corresponds to 
a pair |0^,a)). 

The vectors |(/>y,a) are called the Sturmian basis vectors or simply the Sturmians They satisfy certain 
orthonormality conditions which we shall derive below. We should, however, note that the square-integrable 
Sturmians do not generally constitute a complete set of basis vectors of the Hilbert space Q . There are certain 
potentials Vq, such as the Coulomb potential, that lead to a complete set of square-integrable Sturmians ||]. 

Let us first note that the defining equation (|^) does not determine \(j)i,,a) uniquely. This is reflected in 
the presence of the degeneracy label a. What is uniquely determined by Eq. (|^) is the degeneracy subspace 
Til, spanned by 1), \4>v, 2), • • • , iv)}, where ti, is the degree of degeneracy, i.e., the number of linearly 
independent solutions of Eq. (||) associated with a given (admissible) value of (iv Clearly, we can construct 
an orthonormal basis of Tiu and choose the Sturmian vectors |0^,q;) to be the basis vectors. In other words, 
without loss of generality, we can choose to work with the Sturmians a) satisfying 

((/)^,a|(/)^,7) = , (4) 

where da-y denotes the Kronecker delta function. Clearly, any unitary transformation of "Hi, would lead to a new 
set of Sturmians satisfying (^. Therefore, the condition reduces the freedom in the choice of the Sturmians 
\(l)u,Oi), but does not eliminate it. 

Next, we evaluate the Hermitian adjoint of both sides of Eq. (j^), change {v^ a) to (^, 7), and take the inner 
product of both sides of the resulting equation with 10^, a). This yields 

Pti{(l)i„l\yo\(t>u,a) = {(l)ti,l\{£ - HQ)\(l),,,a) . (5) 

We can compute the right-hand side of this equation using Eq. (|). Substituting the result in (||), we find 

(/3m- 7|Vb|0., a) =0. (6) 

If we define 

7l^o|0.,«) , (7) 
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then we can write Eq. in the form 

(0,„7|V^o|0.,a) =iVJ%. . (8) 

Eq. is the desired orthogonahty property of the Sturmians. We can further simplify Eq. (^), by noting that 
the X matrix N^, formed by A^J" is a Hermitian matrix. This means that we can choose \(j)u,a) in such a 
way that N^, is a diagonal matrix. Making this choice, we have 

K'' = KS^^, (9) 
(0^,7 1^0(2^) I = KS^^So.^, (10) 

where iV^, with a € {1, 2, • • • , £1,}, are eigenvalues of the matrix N,j. Since iV^ is Hermitian, are real. 

In summary, we can choose a set of Sturmian vectors | , a) which are eigenvectors of the matrices . 
Therefore, for each value of {\4>y, 1), • • • , forms an orthonormal eigenbasis of N^, in the degeneracy 

subspace Tiy. However, |0i/,Q!) with different values of v are not orthogonal. Instead, they satisfy a modified 
orthogonality condition, namely (|l^) . 

Now, let us expand a solution \E^a) of the Schrodinger equation (^, in a Sturmian basis corresponding to 
a 'solvable' potential Vqj i-e., seek solutions of the form 

00 

u=0 a=l 

where C" are complex coefficients and v is supposed to take discrete values 0, 1, 2, • • •. Note that if the Sturmians 
\(j),j,a) do not form a complete basis, then Eq. ( [Ill ) yields the eigenvectors that belong to the span of \cf)^,a.). 

The Sturmian Approximation of order N is the approximation in which one neglects all the coefficients C" 
in Eq. ( |ll| ) but those with v belonging to a subset Sn+i of nonnegative integers of order + 1. Alternatively, in 
considering the Sturmian approximation of order A^, one confines the range of the indices (of type) to a fixed 



finite set Sn+i- In this way, the infinite sum X^l^o ' ' ' iii Eq. (11) is replaced by the finite sum X^iyes^ 
We shall abbreviate the latter by Yliv- Sn+i may, in principle, be chosen arbitrarily. We will comment 

on this choice in section III. 

Substituting ( pi] ) in the Schrodinger equation (^ and making use of Eqs. (|) and (|), we find 

^X^C^(ii;-£-l/ + /3,l/o)|0.,a) =0. (12) 



Now, evaluating the inner product of both sides of this equation with 10^,7) and using the orthogonality 
relation (|l^), we obtain 

^ £ [(£; - £)r;; - (i^;; - p^n^.^s,^)] c:^o. (13) 
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Here we have introduced 



:= (0,,,7l0.,a) , (14) 

W^;; := (0p,7l^l0.,a) • (15) 

We can express Eq. (^3|) in a more compact form, if we use a single label for the pair {1^,0). Introducing 
:= [v, a) and M := (/x, 7), we write Eq. (ITsh in the form 



J2[iE~£)TMM~SMM]CM = 0, (16) 
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where 



TmM = T'^v > SmM = WmM - PuNj^SmaT , (17) 

Wmx - w^;", N^^K- (18) 

Eqs. (jl^) form a linear system of homogeneous first order algebraic equations for Cj\f. This system has a 
nontrivial solution provided that the determinant of the matrix of coefficients vanishes, i.e., 

det [{E -£)T~^S] = 0. (f 9) 

Here T and S are matrices with entries Tj^j^ and Sj^j^, respectively. 

Solving Eq. (|l|), we can express E in terms of 5, P^, and the Sturmians \(j)^,a). Now, we recall that 
for a fixed choice of Vq, the coupling constants and the corresponding Sturmians 10^,0;) depend on the 
parameter £. Therefore, Eq. ( |l9| ) yields £' as a function of £. Furthermore, substituting the value of i? = E{£) 
obtained by solving ( pj| ) in ( [T6| ) and solving for the coefficients Cj^, we obtain an expression for the eigenvector 
\E, a) that involves £. As we discuss in the following section, the fact that the Sturmian approximation yields 
the eigenvalues and the eigenvectors of the Hamiltonian as functions of a free parameter seems to have been 
overlooked. This is mainly because there is a choice for the parameter £ that simplifies the calculations. 

It should also be emphasized that Eq. ( p^ is an algebraic equation of order + 1. Therefore, in general, 
it has iV + 1 solutions. This can be understood by noting that in the Sturmian approximation one actually 
approximates the Hilbert space by a finite-dimensional vector space. Consequently, the Hamiltonian is replaced 
with a matrix with a finite number of eigenvalues. 



III. Variational Sturmian Approximation 

In general, the accuracy of the Sturmian approximation depends on the following factors. 

1. Choice of Vq: In practice, Vq must be one of the exactly solvable potentials. Therefore, the available 
choices for Vq are few in number. For the potentials V with bound states, we can choose Vq to be a 
harmonic oscillator potential. For example, for the quartic anharmonic oscillator 

k 

V{x) ^ -x^ + ex^ , (20) 

we shall take 

Vo{x)^^x'. (21) 
V{x) = -Ae-^"'/2 , (22) 

Voix) = ^Xex'' - X . (23) 

This will enable us to compare the results of the Sturmian approximation with those of the perturbation 
theory, for in the limit e — s- 0, V{r) — > Vo(r). Note that multiplying Vq by a positive real number does not 
change the results of the Sturmian approximation. This is simply because we can always absorb such a 
number in the definition of /Jy. 



Similarly, for the Gaussian potential 



we shall take 



2. Choice of the Sturmians included in the sum (16): This is also directly related to the choice of 



the potential Vq. If Vq is obtained from V by a, limiting process as in the case of the potentials (|20| ) and 
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(p2[), then a natural choice for the computation of the n-th energy eigenvalue En and the corresponding 
eigenvectors \En,a) is to include the |0i/,a) with ly equal or close to n. In particular, in the Sturmian 
approximation of order zero, we have 

\En,a) = J2C^\'l^n,a). (24) 

a=l 

3. Choice of the parameter £: The conventional choice |^ for the parameter £ is £ = E. This 
simplifies Eq. considerably. The basic idea pursued in this article is the fact that this simplification 
does not necessarily justify the conventional choice for £. 

It is well-known |^ that the eigenvalue equation (|l|) is equivalent to the variational equation 
In other words, the eigenvalues E are the minima of the expectation value 

_ (mm 

and the eigenvectors are the vectors \E) = \E,a) that minimize (H). This observation suggests that the 
most efficient choice for the parameter £ appearing in the Sturmian approximation is the one that minimizes 
E = E{£). Therefore, the most reliable Sturmian approximation is obtained by choosing £ to be a solution of 

If this equation does not have a solution, then one must cither make another choice for the set Sn+i or proceed 
with a higher order Sturmian approximation. 



IV. Variational Sturmian Approximation Using Harmonic Oscillator 
Sturmians 

Consider a quantum system with the configuration space M, a standard Hamiltonian (^, and a real- valued 
potential V = V{x). Suppose that the system has an infinite number of bound states with nondegenerate 
energy eigenvalues En- Here n e {0, 1, 2, • • •} and Eq stands for the ground state. 

Now, consider the Sturmian basis vectors associated with a harmonic oscillator ^ , 

Fo = K,(x)-^a;2. (27) 

In order to solve Eq. (0) for this choice of Vbj we introduce 



m 
muj,. 



(28) 
(29) 
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K). 4:4'|0)., (31) 
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where |0)i/ is the normaUzed real ground state vector for a harmonic oscillator with mass m and frequency w^. 
That is 

(:,|0).:=(^)'^%-"-V2. (32) 

In view of the similarity of Eq. (|^) with the eigenvalue equation for the potential (iyV^, we can easily deduce 
that 

£ = hLo,{y+]^), (33) 
10.) = (34) 

where v £ {0, 1,2, • • •}. 

We can invert Eqs. (^3|) and ( p^ to express lou and /Si, in terms of £. This yields 

^ 

?i2A:(2zy + 1)2 ^ ^ 



Substituting Eq. (^ in (g9|), we have 

2m£ ao , , 

= = ■ 3/ 

h^{2v + l) 2v + l ' 

Next, we compute the term P^N^ = Pi^{(t>iy\Vo\(f),y) . We can use the properties of the annihilation operator 
Ci/, namely 

a,\£), = VI.\£-l), , al\£), = VeTT\£+l),. , x = {2a,)-^^^{a, + al), (38) 
and the orthonormality of \£)^ to compute 



,{e\x^\i^),^{2a,)-^ {2v + 1)6l. + ^{v + l){v + 2) 5i,,+2 + v(y - 1) (5^,,_2 . (39) 



In view of Eqs. (27) and (|39D, we obtain, after some remarkable simplifications, 

/3.iV. = I . (40) 

Variational Sturmian Approximation of Order Zero 

For the variational Sturmian approximation of order zero, ly — n, and Eq. ( [l9| ) takes the form 

(£;„ - f )T - 5 = , (41) 

where 

T= (0„ !(/.„) = 1, S^W~(3nNn, W = {(bn\V\(bn) ■ (42) 

According to Eqs. (|4l]), (^), and (|40|), the energy eigenvalues i?„ are given by 

= + I . (43) 

Next, we fix the parameter £ using Eq. (p6|). This requires the computation of dW/d£. We first evaluate 
the variation of W, 

6W = ((5(0„|)l/|0„) + (0„|1/(,5|(/.„)) 

= 2(0„|y(5|(/.„)) 

oo 

= 2^ n{n\V\£)n n{mn)n) . (44) 

^=0 



6 



Here we have made use of Eq. (|3^), the fact that the Sturmians and the potential V are real and |£)„ form a 
complete set of basis vectors. 

We can compute n(^|((5|n)„) using the eigenvalue equation 



{Hn + f3nVo)\j)n^£j\j)n 



(45) 



where £j ~ huJn{j + 1/2). Taking the variation of both sides of this equation and computing the inner product 
with \£)n, we find 

nmo\j)n{Spn) 



£j — £i 



for l^j 



Furthermore, using the fact that the eigenfunctions (x|i?)„ are real, we can easily show that 

„(£|(<5|^)„) = 0. 

Eqs. (|^, (|6l) and reduce the computation of 5W to that of 

k 



(46) 



(47) 



(48) 



We have already computed n(^|a;^|n)„ in Eq. (39). Substituting this equation in Eq. ( ^^ and using Eqs. ( p6| ) 
and (^), we find, after some remarkable cancellations. 



5W 



Now, in view of Eqs. 



y/n{n - 1) n{n\V\n - 2)„ - ^{n + l){n + 2) n{n\V\n + 2)„ 



(49) 



and (49), 



dEn 
d£ 



^/n{n - 1) n{n\V\n - 2)„ - ^/{n + l){n + 2) n{n\V\n + 2>„ 



Substituting this equation in Eq. (E6|) yields 



£ = ^J{n + l){n + 2) „(n|F|n + 2)„ - ^n{n-l) n{n\V\n - 2)„ . 

Note that the right-hand side of this equation also involves £. This is because |^)„ depend on £. 
Using Eqs. (^0|) and ( ^ ) we can express the energy eigenvalue i?„ in terms of V . This yields 

En^ n{n\V\n).^ + ]- [y(;rTT)(^rT2) „(n|F|n + 2)„-/;^(;r^„(n|F|n-2>, 



(50) 



For the ground state n = and Eq. (jS^) reduces to 

Eo = o(0|l/|0)o 



V2 



o(0|l^|2)o 



(51) 



(52) 



Note that the vectors appearing in Eqs. ( pl| ) and ( p2| ) are those of Eq. (|3^) with £ being a solution of 
Eq. (1^. 

Eqs. ( pT| ) and ( p2[ ) are of limited importance. In practice, one obtains the energy eigenvalue £"„ by substi- 
tuting the solution of Eq. (|^) in Eq. (||). 

The variational Sturmian approximation of order zero, as outlined above, is a valid approximation scheme, 
if Eq. (|50| ) has a unique positive solution £ . If such a solution does not exist, one may attempt to construct 
higher order variational Sturmian approximations. As we shall see in section V, for all the potentials that we 
have considered, Eq. (|5^) has a unique positive solution. This is very remarkable, for this equation turns out 
to be an algebraic equation of order three for the quartic anharmonic oscillator and the quartic potential, and 
of order four for the Gaussian potential. 
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Variational Sturmian Approximation of Order One 

In the variational sturmian approximation of order one, the number of Sturmians contributing to the eigenvector 
\E) is two. We shall denote them by |0„) and \<j>m)- 

The matrices T and W are Hermitian 2x2 matrices. They can be written in the form 

T=[ ] , W^i " , (53) 

where we have used the fact that the Sturmians are normalized and introduced 

i:=(0™|0„), -.^ {(j)^\V\(j>^) , w:^ {(f)m\V\(j>n) ■ (54) 



Next, we construct the matrix S. In view of Eqs. (|17|) and (|53|), 

^^fv^-PnN^ W* \ ^^^^ 

\ W Vrn- PrnN„i J 

Note that because the Sturmians for the harmonic oscillator are real- valued, t and w are real- valued functions 
of the parameter £. In particular, the matrices T, W and S are real and symmetric. 

Substituting Eqs. (|3|) and ( p5| ) in the Eq. (p£|), making use of Eq. (^, and simphfying the resulting 
expression, we find 

A{E-^r + BiE-^) + C^O, (56) 

where 

t£ 

A:=l-^^ B:=t{t£ + 2w)-{vn+Vm), C := w„t;™ - (y + w)^ (57) 

Note that the coefficients A, i?, and C are functions of £. Eq. ( |56| ) can be easily solved to express E in terms 
of £. The result is 

The next step is to determine £ using the variational principle, i.e., setting dE/d£ — 0. The resulting 
formulas are complicated and we shall not include them here. 
We conclude this section with the following remarks. 

1. As seen from Eq. (|58|), the first order Sturmian approximation leads to a pair of energy eigenvalues. For 
the potentials which are related to Vo via a limiting process, one expects these two eigenvalues to be those 
labelled by m and n. That is, for n < m, 



En = E_, E^^E+. (59) 



2. In the variational Sturmian approximation of order one, there are two variational equations dE±/d£ = 0. 
It is not clear whether these equations lead to a unique minimum for E±{£) with a positive value for £. 
As we shall show in the following sections, for all the specific examples that we have considered each of 
these equations lead to a unique minimum with a positive value for £. Lack of such a solution may be 
interpreted as the failure of the variational Sturmian approximation of order one. 
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Variational Sturmian Approximation of Order Two 



In the variational sturmian approximation of order two, one uses three Sturmians to expand the energy eigen- 
vectors \E). We shall denote these by \(j)ni) where ni G Sz ■= f^2, f^3}■ 
The matrices T, W, and S are given by 



/ 



T = 



S 



1 tl t* 



I 



tl 1 q 

\ t2 t3 1 J 

V3 - PnaNn^ 



W = 



\ 



where we have used the fact that 



V 



i^rii j 



V3 

Wi V2 W3 
y W2 W3 Vi J 

W2 W3 Vi - PniNn^ J 

are normalized and introduced 



(60) 



(61) 



U!l := {(j)n2\V\(j)n-J, W2 ■■= {(f>nAV\<l>n3), := ((?!)„JT^|0„2). 



(62) 
(63) 



Because the harmonic oscillator Sturmian functions are real- valued, te, ve and W£ are real, and T, W, and S 
are real symmetric matrices. 



In view of Eq. (40), we can write the secular equation ( |l9| ) in the form 

AiE-^f + B{E-^f + CiE-^) + D = 0, 



where 



A 



B 



C ;= 



1 - ^ + 2*1^2^3 , 
3 

- \)VI + 2Uit\ - 2{ht2(3 + t3H2 + t2t3^l) , 

i=l 

3 

viV2 + V2V3 + V3V1 + 2(^1^26 + ^366 + ^266) - '^i^t + '^tevi(.£ 



1 



(64) 

(65) 

(66) 

(67) 

(68) 
(69) 



Eq. ( |64D has, in general, three solutions. The desired eigenvalues are the minima of these solutions corresponding 
to positive values of E. Again for the cases where Y is related to Vb by a limiting process the minima of the 
solutions of Eq. (|64|) correspond to En^ , En^ , and En^ . 



V. Applications 

In this section, we apply our general results to compute the energy eigenvalues of a quartic anharmonic oscillator, 
a quartic potential, and a Gaussian potential. 
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The Quartic Anharmonic Oscillator 

Consider the potential 

V{x) 



■ + ex'' 



(70) 



In order to obtain the energy levels of this potential using variational Sturmian approximation of order zero, 
we need to calculate n{n\V\t)n- We first use Eqs. (pq) to compute 



{n\x^\i)r 



(2a„)-2 3(2^2 + 2n + 1)(5^,„ + 4(n + 1) v/(n + + 2) 5^,„+2+ 
2(2n - l)^n{n-l)5u,^2 + \/{n + l)(n + 2){n + 3)(n + 4) 5i^n+A ■ 



\/{n- 3)(n - 2)(n - l)n(5£,„_4 . 
In view of this equation and Eqs. ( |39| ) and (|7C|), 

(2n+l)fc 3(2n2 + 2n+l)e 



1^ : 



4a„ 



4a? 



MV\n + 2)„ = + !)(" + 2) f + 



(71 + l)e 



,.(„iv-i„-2). = v^if;r^(3i; + e^), 

Next, we substitute Eqs. (|7|) and (^ in Eq. (|5^). Using Eq. (|3^), we then obtain 

£^ - pn£ - g„ = , 

where 



(71) 

(72) 
(73) 
(74) 

(75) 



1, 



("+77) , 9« := 



2m2 



(n+ i)2(lln2 + 9n + 4). 



(76) 



It is not difficult to show that Eq. ( |75| ) has a single positive solution [This is true for any positive Pn and q,, 
given by 

^ „ X 1/3 

-1/3 



where 



(|)-'"(1.VT^-.(|)(1) 



i + VT 



27g2 



8n + 4 
11^2 + 9n + 4 



f , 



108?i e 



2.2 



(77) 



(78) 



The right-hand side of Eq. (^^ is manifestly real and positive for r„ < 1. It is not difficult to check that it 
is also real and positive for r„ > 1. In fact, we can express £ in the form 



(79) 
(80) 



where 

(/)„ := tan"^(Vr„ - 1) . 

Note that for r„ < 1, 0„ is imaginary, but cos((/)„/3) is still real and positive. 

Having fixed the parameter £, we can determine the energy eigenvalues En using Eqs. (^) and (^. We 
first use Eqs.(|3^) and (l79|) to compute 




mk 

cos , 
3 V 3 



Then substituting this equation in Eq. ( |7^ ) and using Eq. (p3[), we find 
/3fc, 



(2n+ 1) 



7 + 3tan2(:^^) 



l + tan'(Si) 



(81) 



(82) 
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In particular, the ground state energy is given by 

--(^)\/| 

In Table 1, we list the numerical values obtained using Eq. (^2|) for the first 10 energy levels of a quartic 
anharmonic oscillator with m = 1/2, k — 2, e ~ 1/10 in units where fi = 1. This table also includes the accurate 
numerical values of Ref. pO| , the values obtained using the conventional Sturmian approximation and the zero 
and first order perturbation theory. [The zero and first order perturbation theory yield 

^(0) = h{k/my^^in +1/2), = ^(0) + 3hh\2n^ + 2n + l)/{Amk) , (84) 

respectively.] The relative difference between the results of the variational Sturmian approximation of order zero 
with the highly accurate numer ical resuhs (E*) of Ref. i.e., the quantity \En - E*\/E*, varies between 
3.38 X 10~^ and 2.77 x 10^'^. For the ground state, this number is 1.53 x lO^'^. Even for the lowest lying energy 
levels where perturbation theory yields reliable results, the zero order variational Sturmian approximation 
produces more accurate values than both the zero and first order perturbation theory. As seen from Table 1, 
the variational Sturmian approximation is better than the conventional Sturmian approximation. 

In the remainder of this section we present the results obtained using the first and second order variational 
Sturmian approximation. The numerical results are respectively presented in Tables 2 and 3. 

As we explained in Section IV, in the variational Sturmian approximation of order one one chooses an 
indexing set 52 consisting of two Sturmians to be included in the expansion of the eigenvector \E). One then 
solves the corresponding secular equation (p^, expresses the solutions E± in terms of the parameter £, and 
finds the minima of E±{£). In general, E± are complicated functions of £. However, it turns out that for all 
the cases that we considered E± has a unique minimum corresponding to a positive value of £. 

In order to choose the indexing set 52, we first note that the Sturmian functions {x\(f>n) with even (respectively 
odd) n are even (respectively odd) functions of x. We expect the energy eigenfunctions of the anharmonic 
oscillator ( [70|) to have the same parity structure as the Sturmian functions. This, in particular, suggests that 
in the calculation of Eq we should take ^2 = {0, 2}. 

For a quartic anharmonic oscillator with m = 1/2, fc = 2, e = 1/10, the first order variational Sturmian 
approximation corresponding to 52 — {0,2} yields Eq — = 1.06614 and E2 = -E+ = 5.76117. The value 
obtained for Eq differs from the accurate numerical value by one part in 10**. It is one order of magnitude better 
than the value obtained using the zero order variational Sturmian approximation. The value for E2 is however 
less accurate. One may argue that the choice made for 52 is appropriate only for the ground state. In order to 
compute E2 using the first order variational Sturmian approximation, one may alternatively choose 52 — {2,4}. 
This choice yields E2 = -E- = 5.74558 and E4 = _E+ ~ 9.6637. Again this value for E2 is an order of magnitude 
better than the value obtained using the zero order variational Sturmian approximation, whereas the value for 
E4 is less accurate. One can also try 52 — {0,4}. As expected, this choice yields a less accurate value than the 
choices 52 = {0, 2} and 52 = {2, 4} for both Eq and £'4. 

For the calculation of the first excited state we choose 52 — {1,3}. Then we find Ei ~ — 3.30922 and 
i?3 — E+ — 8.37284. Once again the first order variational Sturmian approximation of order one with the choice 
52 = {1,3} yields a more accurate result for Ei and a less accurate result for E^. 

In general, in the calculation of the energy levels En with n > 2, there are two alternative choices for the 
indexing set 52. In view of the parity properties of the eigenvectors, these are {n, n + 2} and {n — 2, n}. The 
fact that there is no physical reason to distinguish between these two choices suggests that for these levels one 
should consider the second order variational Sturmian approximation with the choice 53 = {n ~ 2,n,n + 2}. 



7 + 3tan2(^) 



l + tan2(^) 



(83) 



11 



Table 3 includes the results of the second order variational Sturmian approximation corresponding to the 
indexing set S3 = {0, 2, 4}. This approximation yields more accurate values for Eq than the zero and first order 
variational Sturmian approximations. However, contrary to our expectation the value obtained for E2 is less 
accurate than the one given by the zero order approximation and the first order approximation with S2 = {2, 4}. 



The Quartic Potential 

Consider the quartic potential 



V{x)=ex^. (85) 



We can easily obtain the energy levels of this potential using the zero order variational Sturmian approximation 
by simply setting A; = in our formulas for the quartic anharmonic oscillator. Substituting A; = in ([76|), we 
can write Eq. (fq) in the form 



1/3 



(n + i)2(lln2 + 9n + 4) 



1/3 



(86) 



In view of Eqs. (g7|, (y), ([Tg), (|86D, and fc = 0, we have 

a„ = h-'{mhe)'/'( ''"J''"J^ ] , (87) 



2n+ 1 

/ 8n + 4 V^'^ fn 2 15 A ^ 
En = n — H n +1] En 

In Table 4, we present the values obtained using Eq. ( |88| ) for the energy levels of a quartic potential with 
m = 1/2 and e = 1 in units where h = 1. This table also includes accurate numerical results E^ and 
the results of the zero and first order WKB approximation given in Refs. ]ll|, |l|. The relative difference 
\En — Ef\/Ef is about 0.04 for the ground state and ranges between 6.4 x 10~^ and 8.7 x 10^'^ for the energy 
levels E2, Ei, Eq, Eg, Eio and Eie . 

Table 5 includes the results of the first and second order variational Sturmian approximation for Eo,E2, and 

E4. 

The Gaussian Potential 

Consider the Gaussian potential 

F(a;) = -Ae-^"'/2 ,. (90) 
In order to apply the results of section IV to this potential, we write V{x) — V{x) — X where 

V{x) = e^'^^^^) . (91) 

Then as e tends to zero, V{x) approaches to the harmonic oscillator potential ( p7|) with k = Ae. 
Clearly, the energy eigenvalues associated with V and V are related by 

En^En-X- (92) 

In the following we use the zero order variational Sturmian approximation to obtain the ground state energy of 
the potential V. The excited energy levels can be obtained similarly. 
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We first note that for the ground state n = 0, and Eq. (50) for the potential V takes the form 



£ = x/2 o(0|y|2)o = V2 o(0|F|2)o = o(0|e-^"'/2|2)o 



U a; e 



x\2)o dx 



(93) 



We can evaluate the right-hand side of (^^ using the well-known expression for the eigenfunctions of the 
harmonic oscillator, namely (^2|) and 

. s 1 /4 

(94) 



(x|2)o=(^)'^'(2aox^-l)e-"-V2. 
Substituting Eqs. (p2|) and (|9^) in Eq. (p3|) and performing the necessary calculations, we find 



where 

Introducing 
we can write Eq. (|9^) in the form 
where 



£{£+pf = AV 



4m 



77 1=1 + -, 

P 

fiv) V'^ - if -r ^0 , 



(95) 

(96) 
(97) 
(98) 
(99) 



It is not difficult to show that for all r > 0, /(?/) has a single minimum at 77 = 3/4. The minimum is 
/(3/4) = -(27/256 + r) < 0. Furthermore, /(O) = /(I) = -r < and lim^^oo /(??) = 00. Therefore, /(??) has 
a single positive root that is greater than 1. This root is given by 

1 



(100) 



where 



IVl-a + b, a:=^C(l + v/T+C^)'/', 
^C(-l + v/TTC^)^/^ C:=^(4r)V3. 



In view of Eq. (^7|) and the fact that 770 > 1, Eq. (|9^) has a single positive solution, namely 

£=pir,^-l). 
Having obtained the parameter £, we next compute 



o(0|l^|0)o = A 1- Vl-'/*' 



(101) 



(102) 



Here we have made use of Eqs. (|3^), (^7\), (|9l|), (p^), and (101). Substituting this equation and Eq. (101) in 
Eq. ( ^3|) and using Eqs. (|9^) and (|9|), we find the ground state energy of the Gaussian potential ( po|) to be 



Eq — Ef) — X 



2^ 



(103) 



In order to reveal the asymptotic behavior of Eq, we investigate the power series expansion of the right-hand 
side of Eq. (|lO^ ). 
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For r ^ 1, i.e. 



-0, 



For r < 1, i.e., \ 



Therefore, for fixed A, 



-A[l 



-1/4 . 



-1/2 



-3/4 



32 



128 



^1+0(^-5/4) 



)[1 



3?- 



and for fixed e, 



lim E{) 

£^0 + 

lim i?o 

e — >oo 



lim i^o 



= -A, 



lim 

e — ^oo 



lim 



4to^A^ 



lim E'n = — oo 



(104) 

(105) 

(106) 
(107) 

(108) 
(109) 



Clearly, the asymptotic behavior of iJg, as given by Eqs. (106) - ( |109| ), agrees with the qualitative analysis of 
the eigenvalue problem for the Gaussian potential. 

It is not difhcult to see that in the limit e — > perturbation theory provides reliable results. Writing the 
Gaussian potential ( pO[ ) in the form 

V = Vo + 5V, (110) 

with Vq given by Eq. (|2^) and performing the standard calculations [Q, we find that the zero and first order 
perturbation theory yield respectively 



E, 



(0) 



E, 



(1) 



= -A(l-r-i/4) 
= -A(l 



^-l/4)-l/2 ^ _ i ^-1/4 ^ 3 ^_i/2 _ ^ ^_3/4 ^ 

2 8 16 



(111) 
(112) 



Here e'^^ is the ground state energy for the Gaussian potential (|9C|) obtained using the ^-th order perturbation 
theory. 



Comparing Eq. (104) with Eqs. (Ill) and ( |112| ), one finds that in the pertubative region where r ^ 1, the 
variational Sturmian approximation of order zero agrees with the results of the perturbation theory. In fact, 
since by construction Eq is the expectation value of the energy of the Sturmian |(/)o), the fact that i?o < E^^^ 
shows that even in the pertubative region the variational Sturmian approximation of order zero is a better 
approximation than the zero order perturbation theory. By the same reasoning, because Eq > E^^ , the first 
order perturbation theory yields a better result. Note however that the wave function obtained in the first order 
perturbation theory is an infinite sum whereas the wave function in the zero order Sturmian approximation is 
given explicitly. 

Another interesting limit is the delta function limit of the potential V where A ~ a^J e/(27r), e — > oo, and 
V(x) — > ~a5{x). Here a is a fixed coupling constant. In this limit r ^ and the ground state energy is given 
by Eq. (105) according to 



Eq 



This result has the same order of magnitude as the exact result: 

-,2 



Eo^ — 



ma 
2n^ 



(113) 



(114) 
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VI. Discussion and Conclusion 

We have outlined a variationally improved Sturmian approximation and applied our results to the harmonic 
oscillator Sturmians. For these Sturmians we could solve the associated variational problem in the zero order 
Sturmian approximation exactly. We have used our variational Sturmian approximation in the calculation of 
the energy levels of various potentials. We have shown that using a few harmonic oscillator Sturmians, one 
obtains quite reliable results. In general, the variational Sturmian approximation is a better approximation 
than the conventional Sturmian approximation. 

Because the harmonic oscillator potential is a confining potential, we expect that the method is more suitable 
for the confining potentials such as the quartic anharmonic oscillator and the quartic potential. We can base 
this argument on a more quantitative reasoning by addressing the problem of classifying the potentials for which 
the Sturmian approximation is exact. It is not difficult to show that these potentials satisfy 



where E, E and C" are constants and :— (x\(^^^ct). Eq. ( |115| ) follows from Eqs. (|^), (|^), (^), and (]ll]). 

For example, the potentials for which the first order harmonic oscillator Sturmian approximation with 
S2 — {0, 2} yields an exact eigcnfunction are of the form 



where ao is a real parameter with the dimension of (length)"^ and ^ is a dimensionless real parameter. As 
seen from Eq. ( |116| ), these potentials tend to the harmonic oscillator potential for |a;| 00. In particular, 
as I a: I ^00, y — > 00. This asymptotic behaviour is also valid for higher order harmonic oscillator Sturmian 
approximations. This observation shows that the harmonic oscillator Sturmian approximation is more reliable 
for confining potentials. 

We conclude this paper with a couple of remarks. 

1. The variational principle used in the variational Sturmian approximation leads to an algebraic (nondif- 
ferential) equation for the parameter £. The acceptable solutions for this equation are those which are 
real and positive. The fact that for all the cases we consider there is a unique real positive solution cor- 
responding to each eigenvalue £"„ is quite remarkable. This observation may be viewed as a consistency 
check for the Sturmian approximation. 

2. In our selection of the Sturmians in the first and higher order Sturmian approximation, we used the 
information about the parity properties of the Sturmians and the energy eigenfunctions. For example we 
ruled out the first order variational Sturmian approximation with S2 — {0,1}. If we perform the necessary 
calculations, we find that for this choice the functions t and w vanish identically and the matrices T and 
S are diagonal. Therefore, the secular equation ( |l9|) yields the same results as the zero order Sturmian 
approximation. This can also be seen from the results of Ref. . 




(115) 




(116) 
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n 


E* 

n 


En 




^CSA 




e'^ 








Et 






El* 





1.065286 


1.06692 


1.5 X 10-3 


1.07500 


9.1 X 10-3 


1.000 


0.061 


1.075 


9.1 X 10-3 


1 


3.306872 


3.31182 


1.5 X 10-3 


3.37500 


0.021 


3.000 


0.032 


3.450 


0.043 


2 


5.747959 


5.75052 


4.5 X 10-4 


5.97500 


0.040 


5.000 


0.13 


5.975 


0.039 


3 


8.352678 


8.34985 


3.4 X 10-4 


8.87500 


0.063 


7.000 


0.16 


8.875 


0.063 


4 


11.09860 


11.0881 


9.5 X 10-4 


12.0750 


0.088 


9.000 


0.19 


12.08 


0.088 


5 


13.96993 


13.9499 


1.4 X 10-3 






11.00 


0.21 


15.58 


0.11 


6 


16.95479 


16.9235 


1.8 X 10^3 






13.00 


0.23 


19.38 


0.14 


7 


20.04386 


19.9998 


2.2 X 10-3 






15.00 


0.25 


23.48 


0.17 


8 


23.22955 


23.1715 


2.5 X 10-3 






17.00 


0.27 


27.88 


0.20 


9 


26.50555 


26.4322 


2.8 X 10-3 






19.00 


0.28 


32.58 


0.23 



Table 1: First 10 energy levels of the Hamiltonian H = + + in units where Ti — 1. E^ are the highly 
accurate numerical values of Ref. ]lO| . £"„ are the values obtained using the zero order variational Sturmian 
approximation. E^^^ are the values obtained by the zero order conventional Sturmian approximation in Ref. 1^. 
Ei°^ and eP are the energy eigenvalues obtained using the zero and first order perturbation theory, respectively. 



52 


Eo 


SEo 


El 


5Ei 


E2 


SE2 


E3 


SE3 


E4 


5E4 


{0,2} 


1.06614 


8.0 X IQ--* 






5.76117 


2.3 X 10-3 










{1,3} 






3.30922 


7.1 X 10-"* 






8.37284 


2.4 X 10-3 






{2,4} 










5.74558 


4.1 X 10-* 






9.66370 


0.13 


{0,4} 


1.06620 


8.6 X 10-* 














9.64502 


0.13 



Table 2: Energy levels of the Hamiltonian H = + x"^ + x^/lO obtained using the first order variational 
Sturmian approximation. 6En stands for \En — E^\/E^. 



n 


E* 


En 


Er, 





1.065286 


1.06613 


7.9 X 10-4 


2 


5.75052 


5.75275 


8.3 X 10-4 


4 


11.09860 


9.68483 


0.127 



Table 3: Energy levels of the Hamiltonian H = + + 2;4/10 obtained using the second order variational 
Sturmian approximation with the choice {0, 2, 4} for the indexing set 1S3. E^ are the highly accurate numerical 
values of Ref. [0. 
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n 


E* 


En 








^(1) 


\l J M n 1 H I TT— 

j^WKB[l) _E^ 


E* 


E* 


E* 





1.060362 


1.10243 


0.040 


0.87 


0.17 


0.98 


0.076 


1 




3.86929 












2 


7.455697 


7.46048 


6.4 X lO"'' 


7.4140 


5.6 X 10-3 


7.4558 


1.4 X 10-^5 


3 




11.6007 












4 


16.261826 


16.1691 


5.7 X 10-3 


16.233615 


1.7 X 10-3 


16.261937 


6.8 X 10-^ 


6 


26.528471 


26.3349 


7.3 X 10-3 


26.506336 


8.3 X 10--* 


26.528513 


1.9 X 10-^ 


8 


37.923001 


37.6218 


7.9 X 10-3 


37.904472 


4.9 X lO-'' 


37.923021 


5.3 X lO-'^ 


10 


50.256255 


49.8404 


8.3 X 10-3 


50.240152 


3.1 X 10-4 


50.256266 


2.2 X lO-'^ 


16 


91.79806 


91.0012 


8.7 X 10-3 











Table 4: Energy levels of the Hamiltonian H — + va units where h = 1. Ef are the highly accurate 
numerical values of Refs. |l^. En are the values obtained using the zero order variational Sturmian approx- 
imation. E^^^^^^ and i?^^^^^^ are the values obtained using the zero and first order WKB approximation 
[ pi] , respectively. 



N 


Sn 


Eo 


\Eo-E*\ 


E2 


\E2-Ef\ 


Ei 


\Ei-Ef\ 


En 


E* 


TP* 

E^ 


2 


{0,2} 


1.08110 


0.0196 


7.60884 


0.0205 






2 


{2,4} 






7.42669 


3.89 X 10-3 


16.4461 


0.0113 


2 


{0,4} 


1.08166 


0.0200 






16.4114 


9.12 X 10-3 


3 


{0,2,4} 


1.08010 


0.0195 


7.56528 


0.0147 


16.5670 


0.0188 



Table 5: Energy levels of the Hamiltonian H = + x"^ obtained using the first and second order variational 
Sturmian approximation. N is the order of the approximation. E^ are the accurate numerical results reported 
in Ref. 0. 
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